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Abstract 

To extend the realm of application of the well known controller design technique 
of interconnection and damping assignment passivity-based control (IDA-PBC) 
of mechanical systems two modifications to the standard method are presented 
in this article. First, similarly to [ij, it is proposed to avoid the splitting of the 
control action into energy-shaping and damping injection terms, but instead to 
carry them out simultaneously. Second, motivated by [^ , we propose to consider 
the inclusion of generalised forces, going beyond the gyroscopic ones used in 
standard IDA-PBC. It is shown that several new controllers for mechanical 
systems designed invoking other (less systematic procedures) that do not satisfy 
the conditions of standard IDA-PBC, actually belong to this new class of SIDA- 
PBC. 

Keywords: Stability of nonlinear systems, passivity-based control, mechanical 
systems. 


1. Introduction 

Stabilization of underactuated mechanical systems shaping their potential 
energy function, and preserving the systems structure, is a simple, robust and 
highly successful technique first introduced in Q. To enlarge its realm of ap¬ 
plication it has been proposed to modify the kinetic energy of the system as 
well. This idea of total energy shaping was first introduced in with the two 
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main approaches being now: the method of controlled Lagrangians Q and In¬ 
terconnection and Damping Assignment Passivity-Based Control (IDA-PBC) 
[1] , see also the closely related work Q . In both cases stabilization (of a desired 
equilibrium) is achieved identifying the class of systems—Lagrangian for the 
first method and Hamiltonian for IDA-PBC—that can possibly be obtained via 
feedback. The conditions under which such a feedback law exists are identified 
by the so-called matching equations, which are a set of quasi-linear partial differ¬ 
ential equations (PDEs), that are naturally split into kinetic energy (KE-PDE) 
and potential energy (PE-PDE). 

Although a lot of research effort has been devoted to the solution of the 
matching equations—see [1, 0 for a recent survey of the existing results—this 
task remains the main stumbling block for the application of these methods. The 
solution of the KE-PDE is simplified by the inclusion of gyroscopic forces in the 
target dynamics, which translates into the presence of a free skew-symmetric 
matrix in the matching equation that reduces the number of PDEs to be solved. 
Due to its Hamiltonian formulation, this term is intrinsic in IDA-PBC, and was 
added to the original controlled Lagrangian method of @,0 —for the first time 


0 —and adopted later ir 


In [ll| it is shown that the PDEs of the 


(^tended) controlled Lagrangian method and IDA-PBC are the same, see also 

m 

Recently, in 3 it has been proposed to consider a more general form for 
these forces, relaxing the skew-symmetry condition. It is claimed in that 
the inclusion of these forces reduces the number of KE-PDEs, but as shown in 
Q this claim turned out to be wrong. One of the objectives of this paper is 
to show that, even though the number of PDEs is not reduced, the inclusion 
of generalised forces effectively extends the realm of application of IDA-PBC. 
A second modification to IDA-PBC proposed in the paper is to simultaneously 
carry out the energy shaping and damping injection steps—instead of doing 
them as separate steps. This modification has been previously reported in [l[, 
where it is shown that the partition into two steps of the design procedure 
induces some loss of generality. In particular, it is shown that (two-step) IDA- 
PBC is not applicable for the induction motor, while SIDA-PBC does apply. 

In the paper we also show that several recent controller designs that do not fit 
in the standard IDA-PBC paradigm, actually belong to this new class of SIDA- 
PBC with generalised forces. In this way, it is shown that these controllers, 
that were derived invoking less systematic procedures, are obtained following 
the well-established SIDA-PBC methodology. 

The remaining of the paper is organized as follows. Section [5] briefly recalls 
the IDA-PBC methodology. Section [3] contains the main result, which is the 
definition of SIDA-PBC with generalised forces. Two recently reported con¬ 
troller design techniques are shown to belong to this class in Section [D The 
paper is wrapped-up with concluding remarks in Section O 


Notation. is the nxn identity matrix and 0„xs is an nx s matrix of zeros, 0„ 
is an n-dimensional column vector of zeros. Given G K, i £ n := {1,..., n}, 
we denote with col(ai) the n-dimensional column vector with elements ai. For 


2 



any matrix A € (A)i G K" denotes the i-th column, {Ay the z-th row 

and (A)ij the zj-th element, G R", z S h, is the Euclidean basis vectors. For 
X G R", S G R”^", 5 = > 0, we denote the Euclidean norm |a;p := x^x, 

and the weighted-norm ||a;||| := x^Sx. Given a function / : R" —>■ R we define 
the differential operators 


V,/ 


dx 





where Xi G R^ is an element of the vector x. For a mapping g : R" —>■ R'", its 
Jacobian matrix is defined as 


Vg:= 


(Vffi)T 

i^gmV 


where gi : R" —>■ R is the z-th element of g. When clear from the context the 
subindex in V will be omitted. To simplify the expressions, the arguments of 
all mappings will be omitted, and will be explicitly written only the first time 
that the mapping is defined. 


2. Standard Interconnection and Damping Assignment PBC 


To make the paper self-contained a brief review of IDA-PBC is presented 
in this section. IDA-PBC was introduced in to control underactuated me¬ 
chanical systems described in port-Hamiltonian (pH) form by 


E : 


q 


. P . 



Jnxn 

In On 


yH{q,p) 


^nxm 

G{q) 


( 1 ) 


where q,p G R" are the generalized position and momenta, respectively, u G R™ 
is the control, G: R" —>■ R"><"* with rank(G) = m < n, the function H: R" x 
R” R, 

H{q,p) ■■=]^P^ M-^{q)p + V{q) (2) 

is the total energy with M: R” —?► R"^”, the positive definite inertia matrix 
and V: R" —>■ R the potential energy. The control objective is to generate 
a state-feedback control that assigns to the closed-loop the stable equilibrium 
(q^p) = (q*y 0), 9* S R". This is achieved in IDA-PBC via a two step procedure. 
The first one, called energy shaping, determines a state-feedback to match the 
pH target dynamics 




q 


OnXn 

M-\q)Md{q) ' 

. P . 



Mq^p) 


\7Hd{q,p) 


with the new total energy function Hd ■ R” x R" —>■ R, 
Hd{q,p) ■■= ^p^ My\q)p + Vdiq), 


( 3 ) 

( 4 ) 
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where Md: K" —>• is positive definite, Vd: K" —K verifies 


g* = arg min 14 ( 9 ), 


( 5 ) 


and J 2 : K." x K” —>■ fulfills the skew-symmetry condition 

MQ:P) = ( 6 ) 


In this case, (( 7 *, 0 ) is a stable equilibrium point of (I3|) with Lyapunov function 
Hd- Indeed, the time derivative of Hd along the trajectories of m takes the 
form 

iJd = pT J 2 V = 0. 

The second step, called -damping injection, is aimed at achieving asymptotic 
stability. This step is carried out feeding back the natural passive output, that 
is, adding to the energy shaping control a term of the form —KpG ^with 
Kp G R”^" positive definite. With this new term we get 


Hd = -\\G^M-^p\\l^ < 0 . 


Asymptotic stability follows if the output G ^p is detectable [13 1. 


To determine the energy-shaping control we equate the right-hand sides of 
Q and (U) to obtain the so-called matching equations 


VgH-Gu = MdM-^yqHd-J2M-^p. (7) 

As shown in Q these equations are equivalent to the solution of the (p-dependent) 
KE-PDE 


{Vqip^M-^p) - Md M-i V,(p^M4V) + 2 J 2 p} = 0., (8) 

the (p-independent) PE-PDE 

G^{yV-MdM-^VVd} = 0s, (9) 

and the (univocally defined) control 

u={G^G)-^G^ [ygH-MdM-^yqHd + J2M-^p] , ( 10 ) 

where G'^: R” —s:=n — m is a full rank left annihilator of G, i.e., 
G'^G = Osxm a-nd rank(G“’‘) = s. 

The success of IDA-PBC relies on the possibility of solving the PDEs (|S]) 
and (|H|). As shown below, the inclusion of generalised forces affects only the KE- 
PDE ©.therefore in the sequel we concentrate our attention on the KE-PDE 
©. In y a more explicit expression for this equation is obtained as follows. 
First, note that to be consistent with ©, whose remaining terms are quadratic 
in p, the free matrix J 2 must be linear in p. Hence, without loss of generality 
we can take J 2 of the form 

n 

=^el M-^pUi{q), (II) 

2=1 
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where : R" —>■ verify Ui{q) = —Uj(q) . To streamline the presentation 

of the result of Q we denote the columns of G-*- as 


G-^(g) 


( 9 )' 

(9). 


where Ufc : R” —>■ R", fc S s := {1,..., s} is given by Vk ■= co\{vki)- Also, we 
introduce the mappings 


Afc : R” R”^”, Bfe : R” R”''”, Tkj : R” ^ R, Wfe : R” ^ R”""”- 


as 


Ak ■■= Md ^j Md, k&s 

n 

Tfcj := y^Vki {MdM~^)ij, k € s, j € n := {I,... ,n} 
Bk ■= Md \ Tfci J Md, kGs 



'vj U{ 


1 _ 

Wk := 


-k 



1 

1 _ 


1 

1 _ 


T 


k G s. 


The proof of the lemma below is given in [8| . 
Lemma 1. The KE-PDE (O is equivalent to the PDEs 


Bk{q) - Ak{q) = Wk{q), k G s. 


( 12 ) 


Note that the left-hand-side of m is a function of the unknown matrix 
Md (and partial derivatives of its components), while the right-hand-side of 
m is independent of the unknown matrix Md (and partial derivatives of its 
components). Hence the number of free elements on the right-hand-side of (E 
entirely determines the number of KE-PDE’s to be solved. It is shown in [8 
that this number equals 

- s (s-I-1) (s-I-2). (13) 

6 

Also, contrary to the claim in [^, the explicit formula m — given in a different 
form also in 
IDA-PBC. 


14| — shows that there is no ansatz for the determination of J 2 in 
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3. Simultaneous IDA PBC with Generalized Forces 


In this paper, motivated by 0 , we investigate the possibility to extending the 
realm of application of IDA-PBC by considering more general external forces. 
In 0 it is proposed to replace the target dynamics in ([3]) by 


q 


. P . 



Onxn M '^{q)Md{q) 

-Md{q)M-^{q) 0„xn 


^Hdiq,p) 


0 

C{q,p) 


(14) 

where C: K" x R” —>■ R.” is a mapping to be defined. Notice that, to ensure 
Hd is a Lyapunov function of the closed-loop—J.e., Hd < 0—the mapping C 
should satisfy 

P^M;['^C < 0. (15) 

Since St and coincide for the particular choice C = J 2 M^^p, it is clear 
that considering these more general forces enlarges the set of desired closed- 
loop dynamics. 

The matching equation now takes the form 


- - V,(p' M-^p) -VV + Gu = -MdM 


-1 


iv,(p^M-ip) + VFd 


(16) 


the KE-PDE ([ 8 ]) becomes 

G^{Vg{p^ M-^p)-MdM-^Vgip^ M-^p) +2 0} = Os, (17) 

while the PE-PDE Q remains unchanged. Stemming from the equation above 
we have two important observations regarding O . 

Ol. Since G(g,0) = 0„ must be satisfied, O can always be expressed in the 
form 

C{q,p) = A{q,p)M;[\q)p, 
for some mapping A : R" x R” —>• 

02 . O must be quadratic in p —this in contrast to the case of J 2 that is linear 
in p. For convenience, and without loss of generality, we take it of the 
form 


with Qi 


2 Oiq,p) = {p^M^ \q) Q,{q) \q)p) a 

i=l 

—>• R"^" free matriees. Consequently, we have 
1 ” 


(18) 


Two consequences of the remarks above are, on one hand, that the target 
dynamics St can be written in the familiar form 


St ■ 


q 


Onxn 

M-^Md ' 

. P . 



A 


ViLd, 


(19) 


6 














and the stability condition pSH now becomes 

M~\q)A{q,p)M~\q)p < 0. (20) 

A sufficient, but not necessary, condition for (12011 to hold is clearly 

A +< 0. 

Notice that, in contrast with the two step design procedure of standard IDA- 
PBC, in this new formulation the energy shaping and the damping injection are 
carried out simultaneously. This is in the spirit of [l[ where it is shown that the 
partition into two steps of the design procedure induces some loss of generality. 
On the other hand, it is easy to see (see Q), that new KE-PDE becomes 

n n 

^ [(i;^ MdM-^ei)V^Md-{vl e,) M^] =-Y^e,vJ Q,(g), 

2=1 2=1 

( 21 ) 

with k G s, and the control law takes the form 

u = {G^ G)-^G^ [VgH - MdM-^VqHd + AM-^p] . ( 22 ) 

Similarly to classical IDA-PBC, the presence of the matrices Qi allows us to 
reduce the number of PDE’s to be solved. Interestingly, this is equal to (fOl) . 
that is, the number of PDEs of IDA-PBC; see Q- In spite of this fact, we show 
in the next section—via a series of examples—that SIDA-PBC with generalised 
forces is applicable to a larger class of systems than standard IDA-PBC. 

We wrap-up this section with a simple proposition that summarises the 
developments presented above and whose proof follows verbatim the proof of 
stability of standard IDA-PBC [^. 

Proposition 1. Consider the underactuated mechanical system ([T]) in closed- 
loop with the control (1221) verifying the following conditions. 

i) Hd and A are given by (|1]) and (fT^ . respectively. 

ii) Md and A satisfy (1^ . 

iii) Aid, Vd and Q verify the matching equations ([H]) and (1^ . 

iv) Md is positive definite and 14 satisfies ©• 

The closed-loop system takes the form m and it has a globally stable equi¬ 
librium at the desired point {q,p) = (g*,0), with Lyapunov function Hd- The 
equilibrium is globally asymptotically stable if 

yo := {A + A^)iM;[^p 

is a detectable output of the closed-loop system. 
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4. Examples of SIDA PBC with Generalised Forces 

In this section we prove that several stabilising controllers for mechanical 
systems—that have been derived invoking other considerations—actually belong 
to the class of SIDA-PBC with generalised forces presented in the previous sec¬ 
tion. More precisely, we prove that replacing the aforementioned state-feedback 
laws in the system © yields the desired target dynamics (HU), j.e., that the 
matching equation (USD holds. 

The definition below is instrumental to articulate our results. 


Definition 1. A state-feedback control law u : K." x K" — >■ R™ for the mechan¬ 
ical system m is said to be a SIDA-PBC with generalised forces if the following 
identity holds true 

-VqH{q,p) + G{q)u{q,p) = -Md{q)M~\q)VqHd{q,p)+A{q,p)Mj^{q)p (23) 

where Hd is of the form o, for some Md positive definite and Vd, A verifying 
(O and (1^ . respectively. Such controllers ensure that the closed-loop system 
takes the pH form (I19II and verify the conditions of Proposition [TJ 


4-1. Energy-shaping without solving PDFs: The controller of 

In 0 a static state-feedback that assigns the Lyapunov function Hd o 
for a class of mechanical systems was given. This control law does not satisfy 
the matching equation o, therefore is not an IDA-PBC. However, we show in 
this subsection that it does satisfy (ESI) — proving that it belongs to the class of 
SIDA-PBC with generalised forces. 

The design of [9| proceeds in two steps, first, a partial feedback linearization 
inner loop is applied to transform the system into Spong’s Normal Form 0. 
Invoking Proposition 7 of [l^ , conditions on M and V are imposed to ensure the 
partially linearized system is still a mechanical system. A consequence of the 
latter is the identification of two new cyclo-passive outputs based upon which 
the controller is designed in a second step. The derivations in 0 are done in the 
Lagrangian form, to fit it into the framework of this paper, we present below its 
pH formulation. 

Consider a mechanical system m with input matrix of the form 

_ An 

Osxm 


Partition the generalised coordinates as g = co\{qa,qu), with € R™ and g^ G 
R'*, which correspond to the actuated and unactuated coordinates, respectively. 
The inertia matrix is conformally partitioned as 


M(g) 


maa{q) m au (9) 1 
mj„(g) m UU (q) J’ 


(24) 


where niaa ■ R" —>■ R"*^"*, ruau '■ R" —>■ and muu '■ R" — 

In Proposition 7 of [l^ it is shown that the mechanical structure is preserved 
after partial feedback linearization if the following conditions are satisfied. 
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Al. The inertia matrix depends only on the unactuated variables qu, *.e., 
M{q) = M{qu). 

A2. The sub-matrix niaa of the inertia matrix is constant. 

A3. The potential energy can be written as V{q) = Va{qa) + Vu[qu)- 
A4. The rows of the matrix mau{qu) satisfy 

^ Quj i^au)k — ^ quki^a,u)j I Vj ^ fc, j,k £ S. (25) 


Under these conditions the system m in closed-loop with the static state- 
feedback control law 

u = UpL{q,p) +v, (26) 

where Mpl : M" x M" —>■ R'" is the partially linearizing feedback given in [l^ . 
see also Section VII of [1], takes the pH form 


4 

P 


0 In 
-In 0 


VH + 


0 

G{qu) 


(27) 


H{q, p) = ^(9«)P + Vuiqu)- 

where 

M{qu) = 

and 

G(gJ := 

Notice that we have dehned a new momenta via 


-rKniQu) 


(28) 


0 

muuiqu) 


p = 


Pa 

Pu 


M{qu)q. 


To complete the controller design the following additional assumptions are 
made in [^. 

A5. The columns of mau{qu) are gradient vector fields, that is, 

'^{mauT = ['^{mauTV, S (29) 


Equivalently, there exists a function Vn ■ R® —>■ R™ such that 

Vn = -mau{qu)qu- (30) 


A6. There exist constants ke,ka,ku G R, Kk,Kj G R™^™, Kk,Kj > 0 such 
that the following holds. 
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(a) det[K{qu)] ^ 0, Vg„ S K®, where K-.W ^ is defined as 


^{Qu) ■ ^u^k^au{Qu^'^uu(.^u')^aui^u)- 

(b) The matrix 


M^\qu) := 


kekalm + klKk X{qu) 
X^iqu) yiqu) 


(32) 


with X{qu) = -kakuKkmau{qu)m~l{qu) and y{qu) = keky,m~^{qu) + 
klrn~^{qu)rnJ^{qu)Kkrnau{qu)rn~^{qu), is positive definite and the func¬ 
tion 

Vd{q) ■■= kekuVuiqu) + ^WKqa + kuVN{qu)\\K,, ( 33 ) 

satisfies condition ([S]), and the minimum is isolated. 

The following proposition is the main stabilization result of Q. 

Proposition 2. Consider the underactuated mechanical system (1771) with mau, 
muu and 14 satisfying Assumptions A5 and A6. The control v : R” x R” —>■ R™ 
given by 


v{q,-p) = -K~ 


k 

ku^k'^au'^uu^Qu^u ^likaQa T ku^N^ ~Tr^k^au^^uu 


^Qu l’^uu'P^]'Pu ^u^k^q['f^au'^uu^u]'^uulP‘ 
KpK {ka^a Pii) ; 


( 34 ) 


with Kp > {) ensures that the closed-loop system has a globally stable equilib¬ 
rium at the desired point (g,p) = (^*,0) with Lyapunov function 

Hd{q,p) = ^P^MJ^p + Vd{q), ( 35 ) 

with and 14 given by (157) and (1551) . respectively. The equilibrium is 
globally asymptotically stable if 

yN •— ^aPa kumiaui^qu^^^uu^^ulPu 


is a detectable output of the closed-loop system. 

Now, we proceed to prove that the control (IMl) is a GSIDA-PBC with gen¬ 
eralised forces. Towards this end, we first notice that the matching equation 
((551) for the system (157() reduces to 

- iv,(p^M-ip) - + Gu = -iMdM-W,(p^M4ip) - MdM-^V,Vd + 

+AM-^p (36) 
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Hence, we must prove that (IMl) verifies (1551) for some A satisfying (1501) . This 
fact is stated in the proposition below whose proof involves a series of long 
computations, therefore, it is given in [Appendix A[ 

Proposition 3. Consider the underactuated mechanical system (1271) with niau, 
niuu and 14 satisfying Assumptions A5 and A6. The control (1551) is a SIDA- 
PBC with generalised forces and 


A(g,p) = -M4V^[M-ip]+M-iV^[A/4ip]-M4iGA-i 


0 : kuKkmaum^u'^JJm^uPu]-2kuKkyq[mau'rn^uPu]'rn^u \ 

G{q)KpG^{q). 


(37) 


Application to the inverted pendulum on a cart 


To illustrate Proposition [3] we consider here the controller for classical cart- 
pendulum example reported in Q . This is a 2-dof system with potential energy 
given by 


V{qu) = mgicos{qu), 


mass matrix 


Me + m m£ cos(g„) 


mi cos{qu) 


and the input matrix is G = col(l, 0), where qa is the position of the car and 
denotes the angle of the pendulum with respect to the up-right vertical position. 
The parameter is the mass of the car, m is the mass of the pendulum and 
i its length. The control objective is to stabilise the up-right vertical position 
of the pendulum. The system satisfies assumptions A1-A4, thus, after using 
a partial-feedack linearising control (j26l) . the dynamics can be written as in 
dUl), with momentum vector p = col(pa,Pu) = col(()a, :^qu), "ninu = mi'^ and 
ruau = micos{qu)- 


In Proposition [5] was used to derive the (locally stabilising) controller 


V 



(38) 


where Ki = 0, /cq = 1 and 


ke + Kk + kuKkm cos^{qu) 



ke + Kk -^^cos{qu) 

. 7- 7. 7^, 


kekumgi cos(g„). 


K{qu) 

Ma^qu) 

Vd{qu) 
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(39) 



















The conditions of Proposition [5] are satisfied if the controller gains verify 
/tg ^ 0, ku ^ 0; k^ -\- K}^ T k^Kj^Tn cos ^ O 5 
for all Qu € (-f, f). 


Some lengthy, but straightforward, calculations show that the control law 
(IITRI) satisfies the matching condition (IM|) with A, derived from (IH71) . given by 


A(g,p) = ]^Md 


0 


kgkuKk 


sin(g„)p„ 


sin(g^)pa 


Md- 


—mi cos(gu) 


Kp [ 1 —micos{qu) ] ■ 


(40) 


4-2. Lyapunov approach for control of underactuated mechanical systems 

Several works have proposed an approach using direct Lyapunov method for 
control design of underactuated mechanicals system (see e.g. [i 3 ,[Il,[ii). In 


the following, we summarise the main idea proposed in these works. 
Consider a mechanical system with dynamics as follows 


q = M ^{q)p 

P = g{q) + fiq,p) + Gu. (4i) 


This dynamics could result from a change of coordinate or a preliminary feed¬ 
back (or change of coordinates) on the mechanical system ([T]) that may not 
preserve neither Lagrangian nor Hamiltonian structure. Notice that the system 
(SB) coincides with the standard mechanical system dB if 

9 {q) + f{q,p) = -^qH{q,p), (42) 


and M is the inertia matrix. 

To proceed with the design, the Lyapunov function candidate 

'Hdiq, p) = ^P^M'f^{q)p -b Vd(g), (43) 

with Md > 0 and g* = argminVd(g) is proposed. The control law is computed 
to ensure that the time derivative of the (H51) along the dynamics (HD) is negative 
semidefinite. That is. 


Hd = p'M^^[giq) + fiq,p) + gu] + -V'[p'M^^p]M-^p + V'VdM-^p 


= pTAf-i 
< 0 


9 (q) + fiq, p) + gu+ -MdM-^Vj -f MdM-^VVd 


(44) 


Now, define the vector 

C{q, P) := 9{q) + f{q, p)+Qu+ [M^"p]p + MdM'^VVd 
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and, recalling observation Ol, rewrite it (without loss of generality) as 


C{q,f)) = A{q,ip)M/p. 

Replacing the equations above in (Hill yields 

pT7W"^A(g,p)7W-^p < 0, (45) 


Proposition 4. The control law obtained via the so-called direct Lyapunov 
approach is a SIDA-PBC with generalized forces. 

Proof The proof follows noting that, from the derivations above, the control 
law should verify 

9 {q)+fiq,V) + Gu = - MdM'^WVd + AM^^p 

= -MdM-^Vgnd + A{q,p)M^^p, (46) 

which coincides with the matching equation (1231) . if we consider a more general 
class of open-loop dynamics for the momenta. This matching equation together 
with the stability condition (H51) shows that the controller is a SIDA-PBC with 
generalized forces, and the closed-loop takes the form (fTiH) . □□□ 


III us- 


Application to the ball and beam system 

We present here the 2-dof example of the ball and beam solved in 
ing the direct Lyapunov method, and show that the resulting controller is a 
SIDA-PBC with generalised forces. The design in 0 first applies a partial- 
feedback linearizing control and a change of coordinate that allows us to write 
the dynamics of the system as follows 


qa 


qu 


' Pa ' 


. P*^ . 



y/2{e+ql) 

0 1 

0 

2{e+q‘i) "I’m 



0 

- sm{qa) 


+ Gu, 


(47) 

(48) 


where qa is the angle of the beam and qu is the position of the ball on the 
beam. The momentum vector is defined as p = M{qu)q, with M.{qu) = 
diag(\/2(e -I- q^), 1), and Q = col(l,0). The control objective is to stabilize 
the equilibrium g* = (0,0). The Lyapunov function candidate has the form 
with 


M/{qu) 


V + gu — y e + 

— \/e + g2 \/2e + qu _ 


and 


Vd{q) = eV2[l - cos(ga)] + y 


qa 


— r= sinh ^ 
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ITj l is as follows 


The controller proposed in 


u 


+ Qu 

•\/e + 92 


sin(ga) + 


— (^S + Kpy/2e + ^2) 


! - (l-u^d CflpQ, C-upU 

v ^ 


pa + Kp^Je + qlpu, 


(49) 


with the functions 

Qupu ^ QuPa 

2yj2e -\- e -\- 2(e + g^) 

qupa 

2\/2e + <72 \/ e + 92 

and parameters iiT, Kp and e positive constants to be chosen. 

We show in |Appendix B| that the controller satisfies the matching equa¬ 
tion (H51) with 


Cu{q,p) ■■= 
Ca{q,p) ■■= 


Hq,p)-=-^Md 

1 

e 


\/«+^ (e-l-92) 

^\/ £ + 92 

<5\/e + 92 ^\/2e -I- 92 . 


V^2e-r<jgi;„p, 

F+9J1 


0 




(50) 


We now verify that the matrix A defined in (1501) satisfies the stability condition 
((051) . For, we notice that the first matrix in ((501) is skew symmetric. Now, 
factoring the term |, the second matrix can be partitioned as 


\/2e + 92 + €^P + 92 


\/2e + 92 \/e + 92 


■ eKp 0 ■ 

\/e + 92 \/2e + 92 . 


. \/e -1- <72 \/2e + 92 . 


0 0 


This matrix is positive definite because S, Kp and e are positive constants and 
the determinant of the hrst right hand matrix equals e. Therefore, the control 
law (1491) is a SIDA-PBC, and the closed-loop dynamics can be written in the 
form (ITOl) . 


5. Conclusions 

An extension to the well known IDA-PBC method for mechanical systems 
has been reported. It essentially consists of two parts: (i) allowing the pres¬ 
ence in the target dynamics of forces, which are more general than the usual 
gyroscopic ones, and (ii) the proposition of simultaneously carrying out the en¬ 
ergy shaping and damping injection steps—instead of doing them as separate 
steps. These two modifications have been previously reported in Q and [l|, 
respectively. 
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It has been shown that several recent controller designs that do not fit in 
the standard IDA-PBC paradigm, actually belong to this new class of SIDA- 
PBC with generalised forces. In this way, it is shown that these controllers, 
that were derived invoking less systematic procedures, are obtained following 
the well-established SIDA-PBC methodology. 
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Appendix A. Proof of Proposition [3] 


We state first the following lemma, whose proof is established using straight¬ 
forward calculations, that will be used below. 

Lemma 2. Given the matrices and G as in (1321) and (1281) . respectively. 
The following relation holds 


M-^G = 


kalm 


K. 


(A.l) 
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The proof of proposition [3] is divided in two parts. First, we verify that (IMll 
satisfies the matching equation (1331) . Second, we prove that A given in (1371) 
satisfies the stability condition (1201) . 

The matching equation (1231) for the system (1271) is equivalent to 


--M-^V,lp'M-^p]-M-^V,V + M-^Gv = --M-^V,lp' M'^p]- 

- M-^VgVa + M-^AM-^p 

Then, the control law should satisy 


[M-'p]p + M,-W,V - [M,-ip]p - 


-M-^V,Vd + M-^AM-^p 

= [M"'P]P + M-W,V - [M-ip]p - M-^V,V, + 

[kuKkmaum-^Vjjm-^PuJ - 2kuKkVg[maum-^p^Jm-^] p^, - 
[M-ip]p + [M-ip]p - M-^GKpG^M-^p 

= M^WgV - M-WgVd - M^^GKpG^M^^p - 

[kuKkmaum-^Vjjm-^Pu] - 2kuKkVg[maum-^p^]m-^] p„ 


-kakuKkUlaum^ 


keku-m^^ + klm^^mJ^Kkmaum. 


1-1 

UU 




Omxs 

'iekum-^ 




h T 

-kum-^ml^ 
yT r™-i. 


^_i ku 


Kiikaqi + kuVN) - M-^GK-^^Kkm^ 


qJ'^uuPulPu + ^GK kuKkVq[maurn^uPu]rn^uPu - GA>G 


ka^m 

-kum-^mj^ 


kuKk-maum^^Vq^V - 


ka^m 

-kum-^rn^^ 


Ki 


{kaqi + kuVN) - '^GK —KkrUaum^u"^[muuPu]Pu + '^GK ^k^ 

ArfeV,[ma«m“.^p„]m“Jp.u - M~^GKpG^AI~^p 


= -MT'^GK- 


ku 


kuKkmaum^u^q^V + Ki(kaqi + ZcuVat) - -^Kkinaum. 


-1 

UU 


^Qu [^uuPu\Pu H“ kuK}^\/q[7Tlau'^uuPu]'^uuP't^ 


- M-^GKpK^ 


(^apa ku’fTT'au’f^uuP'f^^ 
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kuKklTlaum^u'^q^V + Kj^kaqi + kuVN) - 
[^uuPuIPm ~k Vq [WauTOm; Pu]to„^ Pu 

+KpK^{kaPa-kumau'm~^Pu)^, (A.2) 

where we used the definition of A in the second equality, and the relation (lA.lll 
in the sixth equality. The control law (1341) exactly coincides with the term in 
curly brackets hence it satisfies (IA.2I) and, therefore, the matching equation 


M^^Gv= M-^G{-K-^ 


Now, to prove that A given in (1571) satisfies (1701) . we compute first some 
terms of A as follows 

Ai = 

A2 = 

A3 = 

from which we obtain 

-kakuKk [rriaum-^V^{m-lpu) - 2Vq{maum-lpu)m-^] 

{m-^Pu) - 2Vq{mau'm-lpu)m-l] 

(A.3) 

kakuKkmaum~lV]j im~^Pu) 

-kekum-^V^{m-^Pu) - klm-^mJ^Kkmaum-^V^{mZuPu) 

(A.4) 


Ai = 


At = 


Omxm 

O.sxm 


-MJ^GK-^ 




-M-W^ [M-^p] 

7 


M-^yJ[M-^p] 


kuKkruaum^^V[m^uPu\- 
2‘ky^Kk^q ['^au’^uu Pul'^uu 


As = 


-kakuin^^V^ {KkUlaum^uPu) 


-kakum-^Wj {m-^mJ^KkPa) + 

ke ku V J {m-^ p„) + 

{m-^mJ^Kk-maum-^Pu) 

(A.5) 
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Now, we compute (uni) using (13711 uud (IA.3ll - (IA.5l ) 

(1371) = p^[Ai + A2 + A3 — ^GKpG ^^]p 

r n T r n 

^mxm 

. P“ J [ -fcaA:„m“^V^(Arfcma„TO“^p„) 

2kakuKkV^ {maum~^Pu)m~^ 

-kakum-^y^ + 

+klm-lV^ {m-^mJ^KkrnaumZuPu) 

= 2kaku [p^ q{mau'm‘uuP'u-)^uuPu Pa ^9 Pu] 

-klplm-^ [2Vj{maurn~^Pu)Kkmaurn~^Pu - "Vq {m~^rnJ^Kkmaum~^Pu)] - 

-p^M^-^G/^pG^M^-ip 
= -p^M-^GKpG^M-^p, 

which shows that the condition (1201) is satisfied since Kp > 0. 


Pa 

Pa 


- p ' MT^GKpG ' MT^p 


Appendix B. Matching Equation (HOll for the Ball and Beam. 


From the matching equation (1461) . we obtain that the control law should 
satisfy the following equation 

gu = -g{q) - MdM-^VVd - f{q, p) - [M^^p]p + A(g, p)M^^p 


0 

sin(ga) 


-MdM-^VVd- 


0 0 

9«Pa n 

<- 1 .. 


0 0 
0 s 


p- 


--MdM-W^ lM^^p]p + A(q,p)M^^p 


0 

sin(ga) 

0 0 

0 (5 


- 



0 0 ■ 

9gPa n 

L 2(c+9i) ^ J 

-1 

sin{qa) 





(/■uPg 


9gPg 


(<=+9^) 

QgPa n 

(7+S) ^ 


p- 


6 + Kp^y2e+^ —Kpy'e+l^ 
0 S 


vA+9: 


0 




S + Kp \/2e + —Kp y/e + ql 

0 0 


ggPa 


QuPa 


9gPg 


2\/e+9Si/27+ij 2{e+qi) 2y/e+qly/2e+ql 

0 0 
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0 


_ S + Kp^2e + —Kp^ ^ + <& j, 

0 0 J 

which is satisfied by the control law ((1^ . 
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